Finite size effects and 2-string deviations in the 

spin-1 XXZ chains 



Arpad Hegediis 

Research Institute for Particle and Nuclear Physics, 
Hungarian Academy of Sciences, 
H-1525 Budapest 114, P.O.B. 49, Hungary 



Abstract 

In an earlier work [Ij J. Suzuki proposed a set of nonlinear integral equa- 
tions (NLIE) to describe the excited state spectrum of the integrable spin-1 
XXZ chain in its repulsive regime. In this paper we extend his equations for 
the attractive regime of the model, and calculate analytically the conformal 
spectrum of the spin chain. We also discuss the typical root configurations of 
the thermodynamic limit as well as the 2-string deviations of certain excited 
states of the model. Special objects appearing in the NLIE are also treated 
with special care. 



1 Introduction 



Recently spin-1 chains and quantum field theories related to the 19- vertex model [2] 
attract much attention. The spin-1 XXZ chains deducible from the 19-vertex model 
are interesting because on the one hand their Hamiltonian occurs in large QCD 
as the 1-loop anomalous dimension matrix of the single trace operators containing 
the self-dual components of the field strength tensor [3] , and on the other hand the 
determination of correlation functions in these integrable higher spin chains is still an 
area of active research |1] . Furthermore it is known that in the context of light-cone 
approach [5] the inhomogeneous 19-vertex model with alternating inhomogeneities 
provides an integrable lattice regularization for the A/" = 1 supersymmetric sine- 
Gordon model P, [7] . 

In this paper we investigate the finite size effects and 2-string deviations of the 
integrable spin-1 XXZ chain. It is well known that in the thermodynamic limit the 
anti-ferromagnetic ground state of the model is composed of quasi 2-strings, namely 
of pairs of complex roots having imaginary parts close to The deviations of their 
imaginary parts from the values ±| are called 2-string deviations. According to the 
string hypothesis these deviations should be exponentially small in N [8[[9],|10], but 
it turns out that these deviations are much larger; they are of order 1/A^ [TTl m\ . 
This is because the number of 2-strings is not of order one, but A^/2 minus 0(1) such 
that the centers of the outermost 2-strings tend to infinity as tends to infinity. 
This is why these 2-string deviations are not negligible even in the large A^ limit 
calculations of the physical quantities. To treat correctly the technical difficulties 
coming from 2-string deviations one needs to use the so-called nonlinear integral 
equation technique (NLIE). The NLIE technique was originally introduced in [12] 
where with the help of this technique the finite size effects of the ground states of the 
spin- 1/2 and spin-1 XXZ chains were studied calculating analytically their central 
charges and in the spin-1 case the 2-string deviations in the ground state for the first 
time in a unified NLIE approach. Then the NLIE technique was successfully applied 
for describing finite size effects in various integrable spin chains and quantum field 
theories [H [151 [lEl [HI [HI [191 EO] . 

Recently J. Suzuki [1] derived a set of NLIEs (different from the those of [I2]) to 
describe the excited states of the integrable spin-1 XXZ chain in its repulsive regime 
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(i.e. < 7 < |). In this paper we extend his results also for the attractive regime of 
the model (i.e. f < 7 < f )• With the help of our NLIEs, we determine the operator 
content of the conformal field theory describing the thermodynamic limit of the spin 
chain. Furthermore the NLIE enables us to discuss the typical root configurations 
and 2-string deviations of excited states in the large limit. 

The paper is organized as follows: In section 2 some generalities concerning 
Bethe Ansatz solution of the models under consideration are recalled. In section 
3 the classification of Bethe roots is presented. The auxiliary functions and their 
most important properties are listed in section 4. The presentation of the NLIEs can 
be found in section 5. Section 6 contains the counting equations, and section 7 is 
devoted for the discussion of typical root configurations of the thermodynamic limit. 
Section 8 contains the quantitative analysis of 2-string deviations and in section 9 
the calculation of the conformal spectrum of the spin chain is presented. Section 10 
is devoted to special objects. Finally the summary and perspectives for future work 
can be found in section 11. 

2 The definition of the model 

The integrable spin-1 XXZ chain is defined by its Hamiltonian: 

N 

^ = E (^^^ - (^^^)' + ^°^27 {crt - (^- )') - (2 cos 7 - iKa^a^ + a^a^) - 4sin^(5n') , 
1=1 

(1) 

where 

ai = Si- Si+i = + al, a- = S- ■ S-_^^. 

In the rest of the paper we consider N to be even and we impose periodic boundary 
conditions on the spins 5*^+1 = S^, a G {x, y, z}. 

It is well known that the above spin Hamiltonian can be obtained from the 
transfer matrix of the 19- vertex model [U [9l [TO]. Hereafter we consider the 19- 
vertex model with alternating inhomogeneities in order to get access to the light-cone 
Hamiltonian of the J\f = 1 supersymmetric sine-Gordon model as well. 
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Let Vi ~ C''"^^ be the irreducible SU{2) representation with spin li/2 and let 
Rf-'^^\9) the i?-matrices acting on Vi ^Vj satisfying the Yang-Baxter equations 



These i?-matrices can be obtained by fusion [21] from the well known i?-matrix of 
the six- vertex model, and their explicit form can be found in P, [TU] . 

From these i?-matrices one can define a family of transfer matrices with alter- 
nating inhomogeneities Aj = (— 1)*G: 

T,(0, {A,}) = Tr. (i?if - Ai - z^(l + fc)) . . . - - + k))) . (3) 

Due to the Yang-Baxter relation ([2]) the transfer matrices (E]) form a commutative 
family of operators acting on Vh = ^iLi C^: 

[T,ie,{x,}),Tkie\{Xi})]=o. (4) 

Investigating only the case of alternating inhomogeneities we can write for short 
Tk{0, 6) = Tk{6, {(— 1)*6}). Due to integrability and commutativity all these trans- 
fer matrices can be diagonalized by algebraic Bethe Ansatz The eigenvalues 
of the transfer matrices can be characterized by the solutions of the Bethe Ansatz 
equations: 



where 



H9j-t7r) Q(%-Z7r)' ^ l,---,^, (5) 
$m = (smh-(e - Q) smh-(e + Q)Y^\ (6) 

V TT TT / 

M 

g(^) = nsinh^(^-^,), (7) 

J. J. yj- 

J = l 

where 7 is the anisotropy of the model, M is the number of Bethe roots. We recall 
that = — M is the third component of the total spin of the spin chain. The 
eigenvalues of the first two transfer matrices of the fusion hierarchy can be expressed 
by the solutions of the Bethe Ansatz equations as follows: 

T,ie, e) = He - ^.)^^ + + " 



Q{0) ' ' Q{e) 
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The spin Hamiltonian ([T]) is given by the logarithmic derivative of the homogeneous 
transfer matrix T2{6, 0): 

n = -^logUe,0)\e=-q. (10) 
do 2 

On the other hand the regularized finite volume light-cone hamiltonian of the A/" = 1 
supersymmetric sine-Gordon model can be expressed by the inhomogeneous transfer 
matrix T2{9, 6) : 

e^-iii±p)/^ r^T2{±e±t^,e), (11) 

where a = L/N is the lattice constant, and the continuum limit is achieved by taking 
to infinity along with tuning the inhomogeneity parameter as G = In ^ with m 
being the kink mass. 



3 Classification of Bethe roots 

First of all let us parametrize the anisotropy 7 as 

7 = ^. 0<p. (12) 

In this notation the attractive (f < 7 < f ) and the repulsive (0 < 7 < |) regimes 
of the model correspond to the < p < 1 and 1 < p regions respectively. 

All functions entering the Bethe Ansatz equations ([5]) are periodic with respect 
to 7r(j9+2), thus on the complex plane the domain of Bethe roots can be restricted to 
the strip -ZlM < ini0 < niE+Il, The Bethe roots are either real or come in complex 
conjugated pairs, with the exception of self-conjugated roots with Im6'j = 2l^tH). It 
is well known that the ground state of the model is formed by N/2 quasi 2-strings, 
namely pairs of Bethe roots with imaginary parts being close to ±| [HI El [ID]- The 
root configurations of the excitations can be obtained from that of the ground state 
by adding some other new Bethe roots and removing some 2-strings. We classify 
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the Bethe roots as follows: (See figure 1.) 



1. Inner roots: |Ini6'j| < |, 

2. Close roots: | < |Im6'j| < | + niin(l,p) tt, 

3. Wide roots: f + min(l,p) tt < \lmOj\ < 

4. Self-conjugated roots: |Ini^j| = Zli^bH)^ 

In this classification the quasi 2-strings can be either inner- or close-roots. Having 
derived the NLIE governing the finite size effects of the model we will give a more 
precise classification of the Bethe roots. 

It is worth introducing the concept of effective roots as well. To each root 6j with 
I < |Im6'j| < ^^P^'^'^ we associate an "effective root" 9j by the following definition: 

9j^9j-i'^sign{lm9j). (13) 

This transformation is equivalent to the removal of the middle strip of inner roots 
from the fundamental domain of the original Bethe roots. (See figure 2.) 

The classification of effective roots is as follows 
1. Close effective roots: < |Im^j| < min(l,p)7r, 



2. Wide effective roots: min(l,p)7r < |Im6'j| < 



2 ' 



3. Self-conjugated effective roots: |Im6'j| = ZLi^hil^ 

The fundamental domain and the classification of effective roots arc exactly the 
same as those of the Bethe roots in the 6-vertex model, corresponding to the bosonic 
degrees of freedom of the conformal field theory describing the thermodynamic limit 
of the spin chain. 
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Figure 1: Classification of Bethc roots. Crosses stand for the 2-string holes, empty circles denote 
the type I holes, black dots represent the Bethe roots and /x denotes 7rmin(l,p). 

4 Auxiliary functions and their most important 
properties 

In this paper we study the finite size effects of our model with the help of a set 
of nonlinear integral equations (NLIE). In this technique the contribution of the 
N/2 minus 0(1) number of quasi 2-strings is summed up by the integral terms of 
the NLIE, and one has to deal with only a finite number of holes and other type 
of complex roots characterizing the excitation [Ij. This technique is very efficient 
when one is interested in the spectrum of our model in the thermodynamic limit. 
The first step to formulate the NLIE is the definition of proper auxiliary functions. 
Here we use the auxiliary functions introduced in [H [16] . 
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Figure 2: Classification of effective roots. Crosses stand for lioles, effective roots are represented 
by black dots and /i denotes 7rmin(l,p). 



The most natural auxiliary function a{9) is defined by 



a{e) 

1 + a{e) -- 



^{9 -in) Q{e + m) 

1 Q{e) 



Tlie). 



(14) 
(15) 



- in) Q{e + in 

By the help of a{6) the Bethe Ansatz equations (j5]) can be recasted into the form 

a{ej) = -l, j = l,...,M. 
Written in logarithmic form 

i log a{9j) = 2nlj, j = 1 



,M, 



with the Ij quantum numbers being half-integers. From (fTSll one can see that 
i loga(^) takes the value of 2n times a half-integer number also at the positions of 
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the zeroes of Ti{6). Thus taking i log a{6) on the real axis, one can interpret this 
function as the counting function of real Bethe roots, because the integer part of 
i '"^"'■^ M°g"(^) _j_ I provides the sum of the number of those real Bethe roots and real 
zeroes of Ti{6) which lie in the interval [6, 6']. Using this counting function, the real 
zeroes of Ti{6) can be considered as holes in the distribution of real Bethe roots, 
hence we will call them "type I holes" (See figure 1.). The function a{6) is important 
for us for the determination of real zeroes of Ti {6) and for the determination of wide 
roots in the repulsive regime. 

We saw that the previous auxiliary function is the counting function of real roots 
but, since the ground state of the spin-1 chain is formed by 2-strings, we need to 
define an auxiliary function, which can be considered in the thermodynamic limit 
as the counting function of 2-strings. The definition of such an auxiliary function is 
as follows [U HE] : 



+ t^Me + tf) Q{e - tf) '\ 2, 
m = 1 + m = , ' , .3..^^£ Ue). (17) 



We also introduce the complex conjugate of b{9) and B{6) as auxiliary functions 
and we denote them by b{6) and B{6) respectively. 

In the thermodynamic limit the function 4 log b{9) can be considered as the 
counting function of 2-strings. The argument is as follows: from the definitions 
(I16|17l) one can see that 

llog6(0.-^^) =27r/„ IjeZ+^, j = l,...,M, (18) 



- log b{h,) = 2ixh^ , 4 . G Z + ^, (19) 

where hj is a zero of T2{6). In the large limit the 2-strings can be thought to 
be exact, and thus following from ( |T8|T9l) . in this limit considering the function 
jlogb{6) on the real axis one can recognize that the real part of 1 !£sM^^z!£sM^ _)_ x 
provides the sum of the number of real zeros of T2{6) lying in the interval [6, 6'] and 
the number of 2-strings with 6 < KeOj < 6'. In this way real zeroes of the transfer 
matrix eigenvalue T2{6) can be considered as holes in the distribution of 2-strings 
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hence we will call them simply holes. The holes can be depicted in the plane of 
effective roots, so that they lie on the real axis (See figure 2.). 

For the closure of the NLIE we need to define a TBA type of auxiliary function 
as well: 

An additional auxiliary function has to be defined for the determination of wide 
roots in the attractive regime. Its definition is as follows: 

In the context of the NLIE technique the function h{9) turns out to be the most 
convenient to determine the positions of wide roots. It gives the quantization rule 

h{d, - z|) = -1 

for the 9j wide roots lying on the upper half of the complex plane. 

In order to get the proper NLIE one needs to know the analytic properties of 
the main building blocks (i.e. Q{9),Ti{9),T2{9)) of the auxihary functions (11411211) . 
Based on numerical evidence these are as follows [H [16] : 

1. The number of zeroes of Q{9) with imaginary parts being close to ±z| is 
minus 0(1), and it has of order 1 number of zeroes located outside of this region. 

2. Ti{9) is analytic and non zero (ANZ) in the strip Im6' G [— 7r/2,7r/2] apart 
from a finite number of zeroes located exactly on the real axis: 

T,{hf^) = 0, J = 1,...,N^. (22) 

3. Also T2{9) has only a finite number of zeroes in the strip Im^^ G [— 7r/2, 7r/2] and 
they are also distributed along the real axis 

T2{h,)=0, j = 1,...,Nh. (23) 

The knowledge of the analytic properties of these building blocks enables us 
to determine easily the analytic properties of the auxiliary functions (I14ti2ip which 
provide the starting point for the derivation of the NLIE [H [16]. 
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5 The nonlinear integral equations 



In this section we present the NLIE governing the finite size scahng of the inho- 
mogeneous 19-vertex model with alternating inhomogeneities. Previously in [1| the 
NLIE has been derived for the repulsive regime (1 < p), and now we extend these 
results for the attractive regime (0 < p < 1) of the model as well. Apart from 
some technical subtleties, the derivation of the NLIE in the attractive regime is the 
same as in the repulsive one. Skipping the derivation of the NLIE and referring the 
interested reader to references [Il[16j, we simply present the final form of the NLIE: 

\ogb{e) = Ct + tD{e)+ig,{e)+tgbie) + {G*^\nB){e)-{G*^\nB){e) 

+ \im(K+^~' *lnY)(e) (24) 

log?/(^) = Cy + tgyie) + iK+U^\nBm + iK~U^\nBm, 

where In denotes the "fundamental" logarithm function having its branch cut on 
the negative real axis and we introduced the notation for any function / 

The eqs. fl24l) contain three type of convolutions, one of them is the usual one 
containing integration along the real axis 

oo 

{f*9){x)= j dyf{x-y)g{y), 

— CO 

while the other two ones are defined by integrating on the complex plane along the 
integration contours r(t) and T(t) (t G R): 

(/ *r 9){x) = j dz fix - z)g{z), (/ *_ g){x) = j_dz f{x - z)g{z), 

where the curve r{t) is the complex conjugate of r(t). The continuous non-self- 
intersecting contour r(t) has to fulfill the following properties: 

1. Rer(±oo) = ±oo, 

2. < Imr(t) < min(l,p) 7r/2 Vt G M. 
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The first property means that the curve r{t) goes to infinity in both directions 
and the second one is necessary to avoid the poles of the kernels G{x) and K{x). 
Although the form of the eqs. flMl) can describe correctly the finite size effects of 
the spin-1 chain with any choice of the contour r{t) fulfilling the previous properties 
but, since we would like to treat only few a number of roots and holes characterizing 
the excitations, we impose a third condition on r(t) to be fulfilled, namely: 

3. With the exception of a finite number of Bethe roots: |Ini6'j| < 7r/2+Imr(Re^j), 
or equivalently: 

with the exception of a finite number of effective roots: |Ini6'j| < Imr(Re^j). 

This third condition ensures that we need to treat only of order one number of 
roots and holes to characterize the excitations. (See figure 2.) 
The kernel functions G and K of read as 

— oo 

We also introduce the odd primitives of the kernel functions (see appendix C. for 
the choice of branch cuts of Xk{0)), 

e e 
XW^2./..GW. x.W^2./../a.). (26) 



that are important in writing the source terms containing information on the exci- 
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tat ions: 

j=l j=l j=l 

Mc Mw Msc 

i=i j=i i=i 



am = YxK{e-hf), 



i=i 



/ 7r 

9y{0) = \im gy [9 + i- - if] 



V- 

Nh Ms 

j=i j=\ j=i 

Mc Mw Msc 



Y ^^{9 - Cj) - Y XKIli.0 - Wj) - Y XKiiid - w 



c /' 



j=l j=l j=l 

where the second determination of any function: fn{0) is defined as in ^7] 

, //W + /(^-^vrsign(Ime)) Kp 

\f{0)-f{d~ip^sign{lme)) 0<p<l. 



The objects appearing in the source terms of ( |24l) are as follows 



1. Type I holes: {hf}, j = 1, . . . , A^i 



2. Holes: {h^}, j = 1,...,Nh 



3. "Close source objects": {cj} j = l,...,Mc, which are "close effective roots" 
satisfying the condition: Imr(ReCj) < |ImCj| < min(l,p)7r. (See figure 2. and 3.) 



4. Wide effective roots: {wj}, j = 1, . . . , Mw 



5. Self-conjugated effective roots: {wsi }, j = I, . . . , Mgc- 

There are also two types of special objects in the source terms of our equations. It 
is well known that special objects appear when along the integration contour T{t) 
the function B{d) crosses the cut of the logarithm function [T7] . 
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It can be shown that special objects can always be avoided by an appropriate 
choice of the integration contour. On the other hand it can also be shown that for 
any eigenstate of the Hamiltonian the contour r{t) can be chosen in such a way 
that the NLIE must contain at least one special object. Since we want to give 
a formulation of the NLIE being valid for any choice of the non-self-intersecting 
integration contour r(t) we need to take into account the special objects as well. 

In our case the special objects {sj} and {vj} are such points of the contour T{t) 
at which the function In B{6) is exactly on the cut of the logarithm. I.e. at these 
points the value of B{9) is a negative real number 3. 



Im 




Figure 3: Zoom of figure 2. near the real axis. 

The difference between the two types of special objects introduced above is the 
sign of the derivative of the imaginary part of the function log6(^) taken at the 
values of the special objects. For the objects denoted by {sj} this derivative is 
negative, and for the objects {vj} the derivative is positive. In the rest of the paper 
we will call the objects {sj} ordinary special objects and we will call the objects {vj} 

^The sets {sj} and {vj} denote the complex conjugates of the sets {sj} and {vj} respectively. 
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virtual special objects. These special objects get the attribute "virtual", because as 
we will see later in many cases one can get rid of them by an appropriate choice of 
the integration contour r(t). 

The formulation of the aforementioned conditions in the language of the function 
b{9) reads as 

Imlog6(s,) = 27r/,^,, \b{s,)\ > 1, (Imlog6)'(s,) < 0, j = l,...,Ns, (31) 

lm\ogbiv,) = 2nh^, \biv,)\ > 1, (Imlog6)'(t;,) > 0, j = l,...,N^. (32) 

To provide better understanding of these special objects we demonstrate that at spe- 
cial choices of the integration contour r{t) the previous two types of special objects 
correspond to either holes or to the upper parts of the effective root counterparts 
of the quasi 2-strings. For the sake of simplicity hereafter we will call the close 
effective roots coming from quasi 2-strings, effective 2-strings. When the deviations 
of all the quasi 2-strings are less than imin(l,p)| the integration contour r{t) can 
be deformed infinitesimally close to the upper parts of the effective 2-strings and 
the holes. 

Let us single out the the joth 2-string with positive 2-string deviatioij^: 

^jt = ^io ±^ (I + ^io) . Xj^,6j,eR, 0<(5j„<l. 

Then the corresponding effective 2-string takes the form: ± i6jg. In the 

sequel let xq the position of the upper part of the effective 2-string (xq ) or a hole. 
In this case xq fulfills the quantization condition as follows: 

log6(xo) = 27r/oi, JoGZ+1/2. (33) 

Let us assume that the integration contour T(t) runs infinitesimally close to xo, and 
around this point behaves like r{t) ~ Xq + irj with t] being an infinitesimal real 
parameter. Then a Taylor expansion around Xq yields: 

lm\ogb{xo + ir]) = 2% Iq + r] {Re log b)' (xq) + 0{ri'^), (34) 

Re log6(xo + irj) = —rj (Imlog6)'(xo) + 0{rf). (35) 

^Only 2-strings with positive 2-string deviations can appear in the effective plane and so in our 
equations. 
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Thus in the 1] —>■ hmit, the x = Xq + irj point of the contour r(t) can correspond 
to a special type of object (i.e 1 + b{x) is a real negative number) in two cases: 

1. < ?7 and (Imlog6)'(x) < 0, 

2. T] <0 and (Imlog6)'(x) > 0. 

In the first case x is an ordinary special object satisfying ( 13T1) . and in the second 
case X corresponds to a virtual special object defined by (IH^ . 

To summarize: from the previous simple argument it can be seen that the or- 
dinary special objects may appear when the integration contour runs above the 
positions of the upper parts of effective 2-strings or holes , and the virtual spe- 
cial objects may show up when the contour r{t) runs under the upper parts of the 
effective 2-strings. The ordinary special objects correspond to the special objects 
introduced in the 6- vertex model in pT]. 

At this point we should mention that the best choice for the integration contour 
when it runs above all the upper parts of the effective 2-strings and under the upper 
parts of the other close effective roots, but there may be some cases (mostly in 
the inhomogeneous case) when due to large deviations of some quasi 2-strings, the 
effective 2-strings cannot be distinguished uniquely from other close effective roots, 
thus the differentiation between effective 2-strings and other close effective roots 
depends on our choice. 

According to our choice, the integration contour r{t) separates the effective 2- 
strings from those close effective roots which show up in our NLIE as close source 
objects. This definition is arbitrary and depends on the actual choice of the curve 
r(t), but the NLIE ^IM treats this problem consistently, and the final physical 
result, the energy is completely independent of the choice of this separating curve. 

The driving term bulk contribution in the equation for In 6(6*) reads as 




Cb = iT[6b 



She {0,1}, 



. _N Ni 



- Msc + S + n_ mod 2 



2^ 



S 



n_ = 



p + 2 



p + 2 



and 



Cy = inS + niQ{p-l) {Mw + M.J 
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where [...] stands for integer part and 0(x) denotes the Heaviside function. In 
addition to ( |24|) we need two other equations for the determination of type I holes 
and effective wide and self-conjugated roots. For the determination of type I holes 
we need to know a{6) on the real axis: 

-\oga{e) =i6a{e) + {K*^\nB){e)-{K*-\nB){e)-Cy, < |Im^| < |, (36) 

Nh Ms 

Sa{e) = J2^K{0-h,) + J2ixK{e-v,) + XK{e-v,))-^{xK{0-s,) + XK{6 

j=l j=l j=l 

Mc Mw Msc 

- Y.^K{e- c,) - x^niQ - w,) - J2 XKii{e - wg^). 
j=i j=i j=i 

The function necessary to know for the determination of wide and self-conjugated 
effective roots is as follows, for min(l,p) vr < Im6' < ^^^^^^ : 

\oga{e) = tDn{e) + igui{e) + i9Mi{e) + {Gn*^\nB){e)-{Gn*,\nB){e) 
+ {iK-^^)jj*\nY){e), (38) 

where 

d{9) = { ' (39) 

[l/h{9) 0<p<l. 

The source objects appearing in the NLIE are not arbitrary parameters, but they 
have to satisfy certain quantization conditions dictated by the analytic properties 
of the auxiliary functions (1 14112 ip . These are as follows 



For holes: 

^log6(/i,) = 27r4^,, j = 1,...,Nh. (40) 

For ordinary special objects: 

Imlog6(s,) =27r/,^, \Ksj)\ > I, (Imlog6)'(s,) < 0, j = l,...,Ns. 

(41) 

For virtual special objects: 

Imlog6(t^,) = 27r/,^, |K^,)| > 1, (Imlog6)'(t;,) > 1, j = l,...,iV^. 

(42) 
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For close source objects (only for the upper part of the close effective pair) 
1 

i - ■ ■>■ -J 
For wide effective roots 



log6(4) = 27rJT, j = l,...,Mc/2. (43) 



- loga(«;;) = 27rJ^T, J = l,...,M^/2. (44) 

• For self-conjugated effective roots: 

- loga(w;l(^)) = 27r/^T(.), J = l,...,M,,. (45) 

So far we have determined only the upper part of the complex pairs, but the other 
parts can be determined by simple complex conjugation. 

• Finally for type I holes: 

- \oga{hf) = 27r V), j = 1, ...,iVi. (46) 

All the above quantum numbers /«j 's are half integers. A state is then identified by 
a choice of the quantum numbers {Ih^, Icj, ■■■)■ We also mention that the NLIE itself 
can impose constraints on the allowed values of some of these quantum numbers. 
To get the finite size effects of the spin chain one needs to solve f l2^ supplemented 
by the quantization conditions ( l36lH6|) for the two unknown functions 6(r(t)) and 
y{x). The equations ( !36lH6ll also express the fact that the knowledge of b{9) on a 
single curve r{t) and y{x) on the real axis is enough to be able to express all the 
auxiliary functions on the whole complex plane and to determine all the roots of the 
Bethe Ansatz equations 1^. 

Having the solution of the NLIE, the corresponding energy eigenvalue of the spin 
chain can be expressed by it: 

' j=i j=i j=i 

Mc Mw Msc ■ p 

- - E^^^(^- ) + h / ^^^^ 

i=i i=i i=i 

' ^ deK\e)\nB{e)\. 



27r 
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6 Counting equations 



The counting equations are certain sum rules which give the number of holes in terms 
of the numbers of different types of other source objects. These equations play an 
important role in the determination of the constant terms of the NLIE, and in the 
conformal analysis of the spin chain. In our spin-1 chain there are two important 
counting equations: one expresses the number of holes, the other one expresses the 
number of type I holes by the numbers of other source objects. 

Starting from the counting function i log a{6) and applying the argument of De- 
stri and de Vega [TTj one gets the counting equation for type I holes: 

Ni-2N^ = S + Mi + Msc- MR-2n_, (48) 

where Mi denotes the sum of effective roots with the property | < |Imi9j| < ^i^^^tH^ 
Mpt stands for the number of real Bethe roots, and A^^ means the number of the real 
special objects. Real special objects can be either real Bethe roots or type I holes. 
They are called specials because at the positions of real specials the counting function 
of real Bethe roots is no more monotonically increasing: i.e. i ^loga(x) < 0. 

There is also a parity constraint for the number of type I holes. Namely A'^i 
must be even. The proof of this statement is related to the derivation of the second 
counting equation. 

From the definition (fT6l) of b{9) it can be shown that the difference 4 (log6(oo) — 
log6(— oo)) must be equal to 2tt times an integer. Calculating this difference from 
the equations (Ell), one gets that 4(log6(oo) — log6(— cxd)) = A^ivr + 27r x (integer) + 

2w(NH+2NS-2Ns-Mc-2S-2e{p-l){Mw+Mso)+2N+) 

— -. iVlatchmg the two conditions for any 

values of p is possible only if A^^i is even and the expression divided by p in the last 
formula is equal to zero. This leads to the second counting equation: 



Nh + 2N^ - 2Ns = Mc + 2S + 2Q{p- l){Mw + M,,) - 2N+, (49) 



where A^4 



35 
p+2 



p+2 



. In case of the absence of virtual special objects the form 



of this counting equation is very similar to the one of the six- vertex model [T7]. 
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7 Excitations in the large N limit 



In this section using our NLIE ([21]) we discuss qualitatively the typical Bethe root 
configurations of the spin chain (homogeneous case) in the N ^ oo limit. In this 
limit the 2-string deviations are small hence the integration contour r{t) can go 
above the positions of all the effective 2-strings. Moreover the counting function 
I log 6(6') is dominated by the bulk source function D{9) as far as — InA^ < Re6' < 
InA^. In this intermediate regime In B{6) is exponentially small in N, hence the 
integral terms containing In 5(6') can be dropped. The integrals containing \nY{6) 
cannot be dropped, but their contribution is of order one, so their contribution 
is only next to leading order in the large N limit. Their contribution is relevant 
for the determination of 2-string deviations, but irrelevant for the discussion of 
the qualitative root positions of the large N limit. Hereafter we need to discuss 
separately the attractive and repulsive regimes of the model, because the typical 
root configurations other than quasi 2-strings are completely different. 

7.1 Repulsive regime 1 < p 

According to ( l43l) the positions of close source objects are to be determined from the 
function 4 log 6(6'). The quantization condition (H3l) also tells us that the imaginary 
part of ^\ogb{9) at the positions of close source objects must be equal to zero. 
However the driving term D{6) has large imaginary part in the large limit, which 
can be compensated only if the relative positions of close source objects tend to 
certain singularities of source function gb{d)- The singularities of gb{d) comes from 
the function x(^) which has poles at iivr. 

Hence the cancellation of the (of order A^) contribution to the imaginary part of 
J log 6(6') dictates that the relative positions of close source objects must tend to m. 
Consequently, in the large A^ limit, close source objects form pairs with difference vr 
in thier imaginary parts: 



C = Co i ( - + yu) , C =Co-n--/i), Co,/iGM, < < -. 



On the other hand it is known that close source objects appear in complex conjugated 
pairs, thus they must fall either into 2-string configurations on the effective plane 





TT 
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(/i = 0): c^^ = Co ± , or into "effective quartets" with: c 



.±Ti = 



Co ±i (f =F/i), just 



like in the spin-1/2 XXZ chain [27 \ . 



However our NLIE does not give any constraint on the imaginary parts of the 
wide roots, because there is no bulk driving term in fl38l) . This is because Djj{6) = 
in the repulsive regime. Transforming back the effective roots to the language of 
original Bethe roots, one gets that in the large N limit of the spin chain the possible 
root configurations are as follows 

1. Real roots: lm9j = 0, 

2. Inner roots: < |Im6'j| < |, 

3. Quasi 2-strings: IvnOj ~ ±|, 

4. Close roots with: \vn6j ~ ±7r, 

5. Close quartets with: 6'=^Ti ^ j^k^^ ^ ^ < /i < |, 6^0 G M, 

6. Wide roots with arbitrary imaginary parts, 

7. Self- conjugated roots with: Im^^ = ZLilLtHl^ 

This means that in the large N limit, contrary to the prediction of the string 
hypothesis [SI El [ID!, in the spectrum of Bethe roots at best only 1-, 2-, and 3-strings 
and anti-strings accompanied by quartets, wide- and inner-roots without constrained 
imaginary parts can be seen. In this regime there are no constraints on the positions 
of wide- and inner-roots, because of the absence of driving bulk source terms in flHUj) 
and (EE]). 

7.2 Attractive regime < j9 < 1 

Contrary to the repulsive regime in the attractive case the bulk driving term Djj{6) 
of log a{6) is no longer zero imposing further constraint on the imaginary parts of 
wide roots as well. Analyzing the NLIE and following the argumentation of [17] , 




20 



the effective roots fall into the configurations as follows: 



1. Arrays of the first kind are effective root configurations containing close source 
objects as well. These type of arrays consist of effective roots of the form: 



9k^ — 9 ± i {fj, — kpn) , 9k^ — 9 ± i {n — fj, — {k — l)p7i) , A; = 0, . . . , 



2p 



with 9 and < /i being real parameters. At certain special values of /i these arrays 
degenerate. There are two degenerate cases: odd degenerate arrays, which contain 
a self-conjugated effective root at 



9sc^9 + i 



.7r(p + 1) 



9 e 



and accompanying complex pairs at 



2p 



and even degenerate ones, which contain complex-pairs of effective roots at the 
positions 

4 = ^±i!^ii^, k^o,...,' 



2p 



These degenerate arrays contain exactly one pair of close source objects. 



2. Arrays of the second kind contain only wide- and self-conjugated effective roots. 
The odd ones contain a self-conjugated effective root 

.7r(p+ 1) 



9sc = 9 + i- 



9 e 



and wide effective-pairs at 



~^^^^.vr0^j2fc + l)p)^ ;t = 0,...,., 0<.< 



while the even ones contain only wide effective-pairs 

4 = ^±^![il^M, k^O,...,s, 0<s< 



2p 



1, 



2p, 



1. 
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The deviations of the imaginary parts of the previous effective root configurations 
from the formulae hsted above are exponentially small in A^. 

Transforming back the effective roots into original Bethe roots we get the fol- 
lowing possible Bethe root configurations in the attractive regime of the model: 



1. Real roots: ImOj = 0, 

2. Inner roots: < |Im6'j| < |, 

3. Quasi 2-strings: huOj ~ ih|, 



4. Arrays of the first kind with Bethe roots at positions: 



5. Odd degenerate arrays of the first kind: 

6. Even degenerate arrays of the first kind: 

.Tx{2-2kp) 



k = 0,..., 



2p 

k = 0,..., 



2p 



dk = e±i- 



k = 0, 



1 ■ 

2^ 



7. Odd arrays of the second kind: 



2 ' 2 
8. Even arrays of the second kind: 

.7r(2 - 2kp) 



< s < 



2p 



-1, 



9k 



9±i- 



k = 0, 



0< s< 



2p 



Finally we mention that the NLIE does not contain the inner roots as source 
objects, thus it does not impose any constraint on the positions of this type of 
Bethe roots. Knowing the solution of the NLIE their positions can be determined 
from the counting function i log a{6) of ( !36l) by imposing the quantization condition 



i log a( 



27ll 



(50) 



where the positions of inner roots are denoted by 6 



in 
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8 2-string deviations 



In this section we discuss the 2-string deviations for the ground state and for some 
simple excited states of the repulsive regime. The method of the calculation of 
2-string deviations relies on the NLIE technique, and can be extended easily in 
principle for all excited states of the model. 

Previously the leading large N corrections of the 2-string deviations of the ground 
state of the spin-1 XXZ chain were calculated analytically in [II1II2]- Now we extend 
these calculations for excited states as well. Let us introduce the counting function 

z,{e) = ^iogb{e), (51) 

and write the quasi 2-strings as 9j = xj ±i{^ + 6j), where xjs are the string centers 
and SjS are the 2-string deviations in imaginary direction. Then following from the 
definition fll6fl7p of b{6) the quantization conditions hold as follows: 

2ttI- 

Z,{x,+^6,) = ^. (52) 

Let us introduce the functions 

ne) = lim (K+^-' * lnY)(e), (53) 

}Cbie) = -{G *^ \nB){e) - -{G *- \nB){e). (54) 



Then using fl24ti30p ^;,(6') can be written as 



z,{e) = Zo{e) + ^Reic{e). (55) 



where 



Z,m = arctan siuh + + W + M^) + Im^W ^ ^^^^ 

It is important to remark that Zo{6) is always real on the real axis. Analyzing 
the formulae (1331) and (I3U|) . one can recognize that in the large N limit the driving 
contribution for Zf,{6) is in Zq{6), which is real along the real axis and the imaginary 
contribution coming from the term T^Re/C(6') is only of order giving only small 
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correction to Zq{6). From now on we assume that the 2-string deviations are small 
{\6j\ ^ 1) thus we can Taylor expand ( l52l) around the 2-string centers: 

^ = Zo{x,) +zS, Z',{x,) + T^Re/C(x,) + |(Re/C)'(x,) + 0{6'), (57) 

Assuming that SjS are of order and requiring the equality up to this order one 
gets: 

27r/ 



^ = Zoixj), (58) 
_ 1 ReJCjx,) 

The last formula (l59l) provides the analytical formula for the 2-string deviations. 
As it can be seen from fl58p the function Zq{6) yields the quantization rules for the 
2-string centers, thus it can be considered as counting function of 2-string centers. It 
follows that the function Z'q{6) is the density function of 2-string centers, and so the 
2-string deviations are inversely proportional to N times the 2-string density just as 
it was experienced in case of the ground state in [TT], hence our approximation is 
really valid as far as the density is not too small (in the middle of the distribution). 

To get concrete formulae for the 2-string deviations of certain excited states one 
has to solve (IMl) in the large N limit. Then the contribution of the function K,b{d) is 
negligible since it tends to zero exponentially in A^. On the other hand in the large 
A^ limit the contribution of K,{9) is of order 1 providing 2-string deviations. 
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Let us discuss a few examples: 
Ground state 

As it is well known the ground state is formed by only 2-strings and using a r(t) 
contour with appropriately large imaginary part there are no source terms present 
in the NLIE ((241). Analyzing the NLIE in the large limit it turns out that 
y{9) = 1 and }C{9) = |ln2. Substituting these results into the formulae 
and (l59l) one gets the well known formula of [H]: 



In 2 

6i = —— coshxo. (60) 

From this formula one can see that the 2-string deviations can be considered to be 
small as long as \xj\ < In A^. 

2-hole states in the repulsive regime 

From ( l49l) one can see that holes can appear only in pairs thus the excitations with 
smallest energies contain 2-holes. According to fHHl) and fjmi) the two holes can be 
accompanied by a pair of close source objects or by a single self-conjugated root and 
for these configurations the number of type I holes is zero (A^i = 0). Let hi and h2 
the positions of the holes and assume that hi < Then from the large A^ analysis 
of the NLIE one gets for all the possible configurations containing 2-holes that 

y{0) = tanh [ ^ J tanh ^ 



Then the evaluation of 1C{6) p| yields 

ic{e) =tn !.Qi(e- ^^^^ ) - Q^ie - hi) - q^ie -h^)}, (6i) 



where 

1 i 
Q^[e) = -— arctansinh(e) - — Incosh(^), (62) 

Q2ie) = -^^-^\ncosh(^-], (63) 



27r 27r V2. 
and X2(6') denotes x(^) of ( |26l) calculated at p = 2. Knowing }C{9) and the actual 
root configuration accompanying the two holes it is easy to calculate Zq{6) in the 
large A^ limit. Here we do not implement this calculation but rather we concentrate 
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on the calculation of the numerator of fl3Ul) . because the sign of the numerator yields 
the sign of the 2-string deviations. Taking the real part of ( l6T|l along the real axis 
one gets 

RelC(6) = - In cosh ( 9 ^— — ^ | In cosh ( -] In cosh ( ^ | . 

^^2 V 2/2 \ 2 J 2 \ 2 J 

(64) 

Analyzing together with ( l59l) it turns out that the 2-string deviations are neg- 
ative between the positions of the two holes (i.e. 6j < 0, if hi < Xj < /i2), and 
the 2-string deviations are positive outside of this region (i.e. 6j > 0, if /12 < Xj or 
Xj < hi). 

Thus we have shown for all the 2-hole states that the 2-string deviations change 
sign at the positions of holes. Relying on earlier numerical studies of ref. [23] we 
assume that this statement persists for multi hole states as well. To summarize, we 
can say that the 2-string deviations can be positive and negative as well, and the 
2-string deviations change their signs at the positions of the holes. For a numerical 
example see figure 4. 

However we should remark that our derivation and so the above statement con- 
cerning the relation between the positions of holes and the signs of 2-string deviations 
is valid only for the middle of the root distribution, where the roots are densely dis- 
tributed (i.e |Im^| < InA^). Outside of this region nothing definite can be said. For 
example in the presence of ordinary special roots, which are typically located at the 
rarely distributed edges of the root distribution, there are holes at the positions of 
which 2-string deviations do not change sign. This can be naively explained by ( l59l) 
saying that at the positions of ordinary special roots (xj) the density of 2-strings is 
negative (i.e. Z'q^Xj) < 0), thus not only the numerator, but also the denominator 
of ( l59l) changes sign at the positions of certain holes. 

Just to mention an example for this case: in the inhomogeneous case there exist 
a numerical example where Ns = 1, A^^^^ = 2 and the two holes induced by the 
single ordinary special object do not change the sign of the 2-string deviations. (See 
figure 5. in appendix C. of ref. [H].) 
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Figure 4: Deviations of the upper parts of 2-strings in a pure 2-hole state {S — 1) a.t p — 2, 
= 64 with hole quantum numbers //^ = —1/2 and = 12 + 1/2. 

9 Conformal limit 

In this section for even we calculate analytically the conformal spectrum of the 
spin chain. It is well known that the energy of the lowest lying excitations form a 
conformal spectrum in the large N limit, i.e. the energy levels behave as 

where c is the central charge and are the left and right conformal dimensions 
of the underlying conformal field theory. In this thermodynamic limit the energy 
levels can be evaluated without solving the NLIE explicitly [12], [TTl [22]. For the 
analytical calculation of the conformal spectrum we follow the approach described 
in detail in [T^. The positions of the sources for N ^ oo can remain finite (central 
objects), or they can move towards the two infinities as ±lnA^ (left/right movers). 
We introduce the finite parts of their positions 6j by subtracting the divergent 
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contribution: 

We denote the number of right/left moving and central objects by N^'^, Np^, M^'°, . . 
etc. For later convenience we introduce the right /left moving and central spin given 
by 

= 1 (^Np^ + 2N^'^'^ - 2Ar|'° - M^'° - 2 e(p - 1) (mJ'° + M^f)^ . (66) 

According to (H9l) they satisfy the sum rule: 

S+ + S- + 3^^ = S~ N+. 

In the thermodynamic limit the NLIE splits into three separate equations corre- 
sponding to the three asymptotic regions. This is why for all the auxiliary functions 
of the NLIE (11411211) we define the so called kink functions as 

F±(^) = lim F(^±lniV). F G {log6, log?/, logo, . . . }. (67) 

In the large N limit these kink functions satisfy the so called kink equations and the 
energy and momentum can be expressed by them. Performing the above kink limit 
on our NLIE the kink equations take the form 

log6±(^) = C, + tCk±±ie^' + igf'2ie)+igi'2i9) 

+ {G*^^\nB±){9)-{G*^^\nB±){9) + {K-^+'*\nY±){9), (68) 

logl/±(^-2^) =Cy + Cy^ + z~g^^2i9) + {K*^JnB^){9)-iK*-JnB^)i9), (69) 
^ogy±{9) = lim log?/± [9 -i^ + ir]) . (70) 

v^r V 2 / 

- loga± (9) = Cy + Cyi + z^S W + (K \nB±){9) - {K \nB±){9), (71) 
lna±(^) = tCi +t~g!^ii9) + (Gu lni?±)(^^) - {Gjj ln5±)(^^), (72) 

where the functional forms of the source functions of the kink functions (16811721) are 
the same as those of the original equations fl241l39p with the difference that they 
contain only the finite part of the left- and right-moving holes and effective roots 
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respectively. % (See for example [6j). The bulky expressions of the constant terms of 
the equations fl681l72p are listed in appendix B. The finite part of the right and left 
moving objects can be obtained from the kink functions by imposing quantization 
conditions very similarly to fj40l ... l46|) . 

After some manipulations [HI [12], [TTl [22] it turns out that in the thermodynamic 
limit the energy and momentum can be expressed by a sum of dilogarithm functions 
with the 6 — > ±oo limiting values of the kink functions in their argument. One 
group of these limiting values agrees with the limiting values of the original auxiliary 
functions at the infinities, namely: 

6±(±oo) = e^^^-^^ 2 cos(7^), i?±(±oo) = e^^^^^ (73) 

sm (7D j 

l/±(±oo) = r±(±oo) = 4cos(7^)' > 0. (74) 

sm (70) 

The other group of limiting values of the kink functions can be determined from the 
kink equations fl68ti72p . They read as follows 

6±(Too) =6±(too) = 0, (75) 



y±{Too) = i-iyy , 6^ = N+ + 2S^ mod 2. (76) 

Using the dilogarithm sum rule of appendix A, and putting everything together one 
gets that the central charge is 3/2 and the conformal weights take the form 

^^ = ^^t + U^ + ^RY + ^± + J±, (77) 



where 



Q± = S^- (78) 

iV± = ^+ ~ + SS^ -liS^ + N+)S^. (79) 
8 2 

J± = t4± T 2/„± ± 2/,± ± 2/,±T ± 2J^±T ± 4±T T 4(i)± (80) 



2 ^ 2 ^ 2 



scJ ' 



■^Left moving kink functions contain only the left moving (indexed by -) objects, and the right 
moving kink functions contain only the right moving (indexed by +) objects in their source terms. 
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and 

4± = If^± , Is±=Yl ^sf ' = V ' ^^^^ = 5Z , • ■ • etc. 

j=i j=i j=i j=i 

furthermore iV+ = A^+mod2, and Nfj are integers depending on the relative posi- 
tions of complex effective roots and they must be determined case by case. In ( 1771) 
the parameter R is the compactification radius defined by 



R ' ^ 



p + 2 



In order that the formulae fl77ti8ip describe the conformal weights of a c = 3/2 
conformal field theory, the difference — A~ must be independent of R. This 
requirement suffices, because in the conformal limit the momentum eigenvalues of 
the spin chain behaves as P ~ (A"*" — A~)/A^, on the other hand the operator e^^ 
is the operator which translates a state by one lattice site. Thus due to the peri- 
odic boundary conditions after unit translations one gets back the original state, 
i.e. e*^^ = 1. It follows that the eigenvalues of the momentum operator look like 
P = ^ X integer hence the difference A"*" — A~ is independent of the compactifica- 
tion radius R. The fact that A+ — A~ is independent of the compactification radius 
implies some constraints on the possible root configurations in the thermodynamic 
limit. The independence of R can be implemented in two ways: 



1. (5+ = Q- in this case is zero. 

2. (5+ = —Q- in this case = S~ 



To summarize: in the thermodynamic limit for states describing the lowest exci- 
tations only such root configurations appear which carry either zero central spin 
{S'^ = 0) or their left and right moving spin is equal {S~^ = S~). 

^Roughly speaking "lowest excitations" means that these excitations give only of order 
correction to the energy of the ground state. In the repulsive regime this means that N'fj = 0, 
because central holes would give an of order one contribution to the energy. In case of the attractive 
regime also the number of the central wide effective roots is zero when we speak about "lowest 
excitations" in the large N limit. 
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Analyzing the formulae (17711811) one can see that the spin of the state can be 
identified with the winding number of the Gaussian part of the c = 3/2 CFT, and 
the parameters 5^ distinguish the Neveu-Schwartz and Ramond sectors. Depending 
on the state under consideration the sum N± + J± can be either integer or half 
integer, but in the Ramond sector they are always integers. Furthermore the 5*° = 
and = conditions guarantee that 6^ = 6~ ensuring the correct appearance of 
the 1/16 terms in the conformal weights. Moreover further analysis of (1771) shows 
that there is a relation between Q± and the winding number S, namely in the NS 
sector {5y = 0): 

Q±eZ if S e 2Z, (82) 
Q±gZ + ^ if SG2Z + 1, (83) 
while in the Ramond sector {6^ = 1): 

Q± e Z + ^ if S e2Z, (84) 

Q±eZ if 5G2Z + 1. (85) 

Putting together the results of the previous analysis of the conformal weights (1771 - 
ISTl) . one can recognize that they can be interpreted within the framework of c = 3/2 
CFT as conformal weights appearing in the modular invariant partition function of 
Di Francesco et al. [ 23] : 

Z{R) = -TT^ I {XoXl/2 + X1/2X0) 

+ (|Xor + 1x1/21') Y +IXi/i6|' Y > q'^-'-'q^"'"' 

nGZ,mG2Z 2n-me2Z+lj 

where 

Km = 2 (I ± y^) 

are the conformal weights of the Gaussian part of the CFT, ri{q) is the Dedekind 
function, and q = e^'^*'^, r being the modular parameter. 

This analytical result agrees with the earlier conjecture of Alcaraz and Martins 
based on numerical investigation of the Bethe Ansatz equations of the spin-1 XXZ 
chain [21]. 
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To close this section we remark that the resuhs for the conformal weights (1771 - 
ISTl) are analytic in the anisotropy parameter 7 in the entire regime [0, 7r/2] although 
the Bethe ansatz patterns for the excited states look qualitatively different for the 
repulsive and the attractive regimes. The situation is quite similar to that of the 
spin-1/2 chain. There, the same derivation that is successful in the repulsive regime 
is applicable in the attractive regime. The key to this is the large freedom in the 
choice of the integration contour F and the possibility to write down non-linear 
integral equations that are algebraically identical to those for the ground state but 
using deformed contours. In this approach there are no additional "driving terms" 
in the NLIEs and the conformal charges are given in terms of dilogarithm functions 
with "non-standard" integration contours [l3]. The analyticity of conformal weights 
in the spectral parameter indicates that this approach can probably be extended for 
higher spin cases too. 



10 Ordinary and virtual special objects 

In this section we would like to shed more light on the appearance of virtual special 
objects in the NLIE, therefore we will consider the ground state as simplest example. 
The ground state of the model is formed by N/2 quasi 2-strings with quite small 
deviations from ±f in their imaginary parts. First let us write down the NLIE when 
the contour is chosen to be a straight line running above all the effective 2-strings 
(i.e. Fi(t) = t + irji with r/i > 6j Vj). (See figure 5.) In this case the NLIE takes 
the form of (Elj) with trivial source terms: 

gii9)=gti9)=gyi9) = 0. 

In this case there are neither ordinary nor virtual special objects present. 

Next we consider another straight contour F2(t) running under exactly two ef- 
fective 2-strings: (i.e. F2(t) = t + irj2 with rj2 < 5i = 6n/2 and ri2 > Sj for all other 
values of j). (See figure 5.) In this case the non-bulk source terms of the NLIE are 
not zero anymore! Two pairs of close source objects corresponding to the two effec- 
tive 2-strings being above T2(t) appear as source terms and also two virtual special 
objects show up in the source term of the NLIE. In the language of the counting 
equation (l49l) the appearance of the four close source objects (i.e Mc = 4) entails 
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Figure 5: Effective plane and graphical demonstration of the free choice of the integration contour 
in case of the ground state of the model. The crosses lying on the contour r2 denote the positions 
of the virtual special objects. 

the appearance of two virtual special objects ensuring that we are still speaking 
about the ground state (i.e. Nh = My/ = Msc = S = 0). Then compared to the 
previous case the source terms of (|2^ are modified as 

9biO) = x{0-vi)+xid-vi)+xid-VN/2)+xid-VN/2) (86) 

- Xid - Cl) - Xid - Cl) - Xid - Cn/2) - Xid - Cn/2), 
~9y{0) = XKiO - Vi) + Xk{0 - Vi) + Xk{0 - VN/2) + Xk{9 - VN/2) (87) 

- XKiO - Cl) - XKiO - Cl) - XKiO - Cn/2) - Xk{,0 - CN/2), 

and gi{0) = 0. As for the notation vi and vn/2 denotes the positions of the vir- 
tual special objects and the set {ci, CAr/2, ci, CAr/2} = {±xi ±iSi} denotes the set of 
close source objects. These source objects are subjected to quantization conditions, 
namely ci, and cn/2 satisfy a condition like (H3l) . and vi, and vn/2 have to satisfy 
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032]) • The other source objects are simply complex conjugates of these ones. This 
was the simplest example how the virtual special objects enter the NLIE. Neverthe- 
less we have to notice that although the two different formulations of the NLIE for 
the ground state are equivalent but they are not equally practical! The second for- 
mulation containing virtual special objects is not appropriate for numerical studies, 
because the usual large iteration of the NLIE fails to converge, whilst the first 
formulation of the NLIE can be solved easily numerically. Certainly when analytical 
calculations are considered (e.g. calculation of conformal weights) they can be used 
equally well, and they give the same result. To summarize: the same state can be 
formulated in many ways by means of NLIE by the modification of the contour r(t). 
The different formulations differ from each other in the number of virtual special 
objects occuring in the source terms. 

The ordinary special objects defined by (1311) according to the counting equation 
f HUj) generally accompanied by the appearance of two holes, just like in the 6- vertex 
model. Taking a glance at the NLIE (12^ . ordinary special objects may pop up in 
those regimes where the derivative of the bulk source term D{6) becomes small. This 
means that in the spin chain they can appear at the edges of the root distribution 
(i.e. InA^ < |Re^|). In case of alternating inhomogeneities if the inhomogeneity 
parameter B scales as In as ^ oo then ordinary special objects may emerge in 
the middle of the root distribution or at its edges. In the inhomogeneous case a nice 
numerical example for ordinary special objects can be found in appendix C of ref. 
[6]. We remark that when any type of special objects show up in the form of the 
NLIE then its numerical solution fails to converge, thus in the presence of specials 
one can rely on only analytical manipulations of the NLIE. 

11 Summary and perspectives 

In this report we extended J. Suzuki's NLIE for the attractive regime (f < 7 < f ) 
of the inhomogeneous 19-vertex model with alternating inhomogeneities. Then we 
analyzed the NLIE corresponding to the homogeneous case to describe the finite 
size spectrum of the integrable spin-1 XXZ chain. Using the NLIE we determined 
analytically the conformal spectrum emerging in the large N limit of the spin chain. 
Our result agreed with the one proposed earlier by Alcaraz and Martins based on 
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numerical investigation of the Bethe Ansatz equations [24j. By means of the NLIE 
we discussed the typical Bethe root configurations of the thermodynamic limit, and 
we worked out a method to determine the 2-string deviations for any eigenstate of 
the model. Moreover we have shown analytically for a class of excited states that, in 
accordance with earlier numerical observations [23j , the holes in the sea of 2-strings 
are the positions were the 2-string deviations change sign. 

Finally, we would like to make a few comments on the possible further applica- 
tions of the NLIE. First of all it is worth mentioning that recently, in the framework 
of quantum inverse scattering method, remarkable progress has been made in the 
calculation of form factors and correlation functions of integrable spin chains [U [26] . 
The correlation functions are represented as complicated multiple integrals, and the 
form factors of local operators are expressed as determinants of matrices, whose ma- 
trix elements are known analytic functions of Bethe roots of the eigenstates involved. 
The multiple integral representations of correlation functions are too complicated for 
their explicit analytical or numerical evaluation. However as it was done in [26], one 
can try to calculate the correlation functions numerically by formulating them as se- 
ries of appropriate form factors. In this case the form factors are determined by their 
determinant representations involving the Bethe roots of the included eigenstates. 

This method requires the knowledge of the Bethe roots for many eigenstates. 
Nevertheless the solution of the Bethe Ansatz equations in higher spin XXZ chains 
is not easy for large because of the string deviations. At this point the NLIE can be 
useful since it can be numerically solved easily for large A^ and using the definitions 
of the auxiliary functions, all the positions of the Bethe roots can be extracted from 
the numerical solution, providing a good starting point for the numerical calculation 
of form factors. 

Another area of applications of the NLIE is the description of finite size effects in 
1+1 dimensional integrable quantum field theories. It is known that the NLIE (12^ 
with appropriately tuned inhomogeneity parameter describe the finite size effects in 
the A/" = 1 supersymmetric sine-Gordon model [6], [7]. The NLIE governing the finite 
size effects of the A/" = 1 supersymmetric sine-Gordon model in the repulsive regime 
was investigated in and the ground state of the model with Dirichlet boundary 
conditions was studied in [7J . The extension of the NLIE technique to the attractive 
regime, presented in this paper, enables one to investigate the attractive regime of 
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the model as well. 
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Appendix A 



In section 9, in the calculation of the conformal weights, the following dilogarithmic 
sum must be calculated: 

sin (375) 



^0 = 2 <! L+ [e^^^^2 cos(75)] + L+ cos(75)] + 

where 



sin (75) J J 



2 y I y 1 + y 



Using the dilogarithm identity of appendix B of [6], it can be easily proven that 
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where 

N+ = N+ mod 2 
and we made the choice of ln(— 1) = zTr. 



TiH - I N+TT In (4 cos2(75)) , (90) 
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Appendix B 

This appendix is devoted to list the constants emerging in the kink equations (16811720 . 
They are as follows 

a± = ±2x00(5 -S+-N+)±n(j- (M,, - M±) + ^ ' \ + 2nK^, (92) 



Cy+ = ivr (5 + Ml - - 5+ - M+) , 



(93) 
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Cy^ = in {S- + Aff + 2M- - S - Mi + N+- 2lVQ , (94) 

C-a± = ±27r {S-S^- N+) + 2nNa±, (95) 

where the integers and Na± depend on the relative positions of the wide effective 
roots, and they must be determined case by case. 



Appendix C 

In this appendix the most important properties of the function Xk{(^) will be clari- 
fied. The function Xk{0) is defined as the odd primitive of the kernel 2t: K{9) and 
given by the formula: 

sinh ii- + 

^""^^^^'^""^M^^' Xi^W = -Xi^M) V^eC. (96) 

The branch cuts are chosen to run parallel to the real axis so that XK{d) be an odd 
real analytic function on the entire complex plane and continuous along the real 
axis. In this case Xk{(^) is not periodic anymore with respect to 27ri. It is periodic 
only modulo 27r, i.e. the following identity holds: 

Xi^(^ + 27ri)=Xi^W-27r. 

It follows that the distance between the consecutive cuts is 2711 and the jump of 
Xk{(^) is equal to —2ti at each branch cut crossed from below to up. The choice of 
branch cuts is depicted in figure 6. 
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Figure 6: Locations of the branch cuts of 
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